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We  report  the  use  of  Yo.7Cao.3Ba2Cu307_5  as  an  epitaxial  barrier  between 
YBa2Cu3C>7-5  (YBCO)  electrodes  in  high-Tc  SNS  edge  junctions.  Ca-doped  YBCO  is 
an  overdoped  version  of  YBCO,  and  it  has  an  excellent  lattice  and  thermal 
expansion  match  with  YBCO.  We  show  that  these  junctions  exhibit  clean  interfaces 
with  resistances  smaller  than  10"10  O  cm2.  We  present  the  temperature  dependence 
of  the  critical  current  density  }c  and  the  junction  resistance  Rn  for  junctions  with  200 
A,  400  A,  and  600  A  thick  barriers  of  Yo.7Cao.3Ba2Cu307-5.  As  the  temperature 
decreases,  the  resistance  of  the  junction  decreases  faster  than  the  resistance 
calculated  from  the  resistivity  of  a  Yo.7Cao.3Ba2Cu307-5  film.  We  analyze  our  data  by 
considering  various  possibilities  for  the  microstructure  in  the  barrier  region.  We 
discuss  the  ideal  conventional  proximity  effect,  diffusion  between  the  barriers  and 
the  electrodes,  and  doping  fluctuations  inside  the  barrier.  We  find  that 
combinations  of  these  ideas  explain  the  observed  behavior  with  reasonable 
parameters. 
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I.  Introduction 


High  temperature  superconductor-normal-superconductor  (SNS)  Josephson 
junctions  have  attracted  considerable  attention  because  of  their  possible  electronic 
applications,  such  as  ultra  fast  digital  logic  circuits  and  shift  registers  or  their  use  in 
SQUIDs.  Josephson  junctions  also  provide  an  elegant  method  to  test  phenomena 
like  the  proximity  effect  [1]  or  to  investigate  such  issues  as  the  wave  function 
symmetry  of  YBCO  [2,  3,  4].  In  order  to  use  these  junctions  as  a  model  system, 
especially  for  investigation  of  a  possible  high-Tc  proximity  effect,  the  junction 
interfaces  must  be  clean,  i.e.  the  interface  resistances  at  the  S-N  boundaries  should 
be  very  small.  The  important  problem  of  interface  resistance  in  high-Tc  junctions 
has  indeed  been  recognized  for  some  time.  For  instance,  junctions  made  with  Au  or 
Ag  [5],  CaRuC>3,  [6]  and  SrRu03  [7]  as  a  barrier  are  known  to  have  interface  resistances 
of  the  order  of  10'8-10'7  H-cm2,  similar  to.  the  values  observed  in  grain-boundary 
weak  links.  These  high  values  influence  the  junction  parameters,  such  as  the 
magnitude  and  reproducibility  of  the  critical  current  IC/  which  limit  the  applicability 
of  such  junctions  in  electronic  circuits.  The  origin  of  the  high  resistance  in  both 
SNS  and  grain-boundary  junctions  is  believed  to  be  oxygen  deficient  or  oxygen 
disordered  YBCO  adjacent  to  the  interface  [8,9].  We  have  shown  recently  that  the 
mismatch  of  thermal  expansion  coefficients  at  the  S-N  interface  may  be  a 
contributing  factor  to  the  excess  interface  resistance  [8].  An  evident  way  to  avoid  this 
problem  is  to  use  a  barrier  material  which  is  both  epitaxial  and  has  a  matching 
coefficient  of  expansion.  This  goal  can  be  achieved  by  using  barriers  which  have  the 
same  layered  structure  as  YBCO  but  with  a  lower  Tc.  The  doped  YBCO  system  is 
thus  an  excellent  candidate. 

In  recent  articles  [8,  10]  we  have  reported  that  junctions  with  barriers  of 
YBa2(Cui-xCox)307-5  and  Y0.7Ca0.3Ba2Cu3O7-5  have  interface  resistances  as  low  as 
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10'10  £Tcm2,  i.e.  two  orders  of  magnitude  lower  than  grain  boundaries  or  junctions 
made  with  Au,  Ag,  and  CaRuC>3  [5,  6].  Those  low  interface  resistances  have  been 
reproduced  in  the  edge  junction  geometry  with  Co-doped  YBCO  by  Koren  et  al.  [11] 
and  in  the  Ca-doped  YBCO  system  by  Jia  et  al.  [12].  Our  previous  study  showed  that 
for  Co-doped  YBCO  junctions,  the  temperature  dependencies  of  Rn  and  lc  are 
different  from  junctions  which  have  higher  interface  resistances  [5].  Furthermore, 
we  have  demonstrated  that  the  IC(T)  curves  can  be  well  described  by  the  framework 
of  de  Gennes'  proximity  effect  theory  [13]. 

It  is  important  to  recognize,  however,  that  these  results  have  only  been 
obtained  with  superconductor  barrier  materials  operating  above  their  Tc. 
Furthermore,  the  Co-doped  YBCO  junctions  utilized  an  underdoped  version  of 
YBCO  as  the  barrier,  and  it  will  be  useful  for  the  understanding  of  SNS  junctions  to 
know  how  the  properties  of  the  junctions  can  be  affected  by  the  carrier  density  of  the 
barrier.  In  order  to  address  this  latter  point,  we  present  the  properties  of  SNS 
junctions  with  barriers  of  Ca-doped  YBCO,  which  is  an  overdoped  version  of  YBCO. 
Although  we  will  show  that  some  of  the  properties  of  these  Ca-doped  YBCO 
junctions  can  be  understood  within  the  de  Gennes  theory,  they  can  also  be  explained 
by  models  invoking  diffusion  and  percolation.  This  follows  from  the  fact  that  the 
barrier  itself  is  a  superconductor  with  a  reduced  Tc,  together  with  the  recognition  of 
a  possible  inhomogeneous  distribution  of  oxygen  and/ or  other  defects  on  a  length 
scale  bigger  than  the  superconductor  coherence  length  a  point  for  which  we 
present  some  evidence. 

II.  Experimental  Procedure 

The  effect  of  doping  YBCO,  especially  with  Ca,  has  been  studied  extensively 
[14].  The  substitution  of  the  Y3*  ion  by  Ca2+  is  expected  to  increase  the  carrier  density 
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by  introducing  one  additional  hole  per  substituted  ion,  provided  that  the 
substitution  does  not  create  compensating  defects  such  as  oxygen  vacancies.  The 
critical  temperature  decreases  with  increasing  doping  and  can  be  lowered  in  thin 
films  to  around  50  K  by  40%  Ca  substitution  [14].  Fig.  1  shows  the  temperature 
dependence  of  the  resistivity  of  an  epitaxial  c-axis-oriented  Yo.7Cao.3Ba2Cu307-5  film. 
This  film  was  deposited  by  laser  ablation  on  a  LaAlC>3  substrate  at  780  °C  in  400 
mTorr  of  O2  and  was  then  cooled  in  600  Torr  of  O2.  X-ray  diffraction  shows  no 
impurity-phase  peaks,  and  the  measured  lattice  constant  in  the  growth  direction  is 
about  11.8  A.  The  resistivity  of  the  Ca-doped  YBCO  is  lower  than  YBCO  by  roughly 
the  amount  expected  from  the  increased  carrier  density;  TC(R  =  0)  is  about  50  K. 
However,  the  zero-resistance  Tc  of  our  unpattemed  and  identically  deposited 
Yo.7Cao.3Ba2Cu307-5  films  varies  randomly  between  40  K  and  60  K.  In  addition, 
Appelboom  et  al.  [14]  have  reported  a  Tc  of  60  K  for  fully-oxygenated 
Yo.7Cao.3Ba2Cu307-s  films  fabricated  by  co-evaporation,  and  Kucera  et  al.  [15]  found 
a  Tc  of  83  K  for  annealed  sputtered  Yo.7Cao.3Ba2Cu3C>7-s  films.  Since  Ca-doped  YBCO 
is  known  to  be  oxygen  deficient,  this  wide  variation  probably  reflects  the  high 
sensitivity  of  Tc  to  the  oxygen  concentration  in  this  material  [15,  16]  .  The  Tc  of  Ca- 
doped  YBCO  films  is  also  quite  sensitive  to  the  oxygen  pressure  used  during  the 
cooling  of  the  sample.  For  example,  cooling  our  laser-ablated  films  in  an  oxygen 
pressure  of  5  Torr  increases  the  Tc  to  ~80  K. 

Fig.  2  shows  the  geometry  of  our  SNS  edge  junctions.  The  YBCO/ Ca-doped 
YBCO/YBCO  junctions  were  fabricated  by  the  same  process  used  for  our  previously- 
described  CaRu03  junctions  [6].  Epitaxial  layers  were  grown  by  laser  ablation  on  a 
YSZ  substrate  with  a  Ce02  buffer  layer.  Thicknesses  of  about  2000  A  for  both  the 
YBCO  electrodes  and  the  SrTiOs  insulating  layer  were  obtained.  Patterning  was 
accomplished  by  Ar  ion-beam  milling  to  obtain  five  4-|im-wide  junctions  per  chip. 
According  to  cross-sectional  scanning  electron  microscope  studies,  the  angle  of  the 
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slope  of  the  first  YBCO  edge  is  about  20°  to  the  plane  of  the  substrate.  In  order  to  test 
the  YBCO  surface  after  the  ion  milling,  we  made  an  edge  junction  structure  with  no 
barrier  layer.  The  critical  current  density  was  higher  than  1  x  10^  A/ cm2  at  77  K, 
suggesting  that  the  ion-beam  damaged  YBCO  surface  was  restored  during  the 
heating  for  the  second  YBCO  deposition.  The  quality  of  the  SrTiOs  insulating  layer 
was  tested  by  making  a  few  junctions  with  a  100  A  thick  barrier  of  SrTiOs-  These 
junctions  did  not  show  any  zero-voltage  current,  even  though  the  resistance  of  the 
junctions  varied  between  tens  of  ohms  and  hundreds  of  kilohms. 

III.  Results 

Fig.  3  shows  the  current-voltage  (I-V)  characteristics  of  a  4-|im-wide  junction 
with  a  400  A  thick  barrier  of  Ca-doped  YBCO  at  several  temperatures.  At  high 
temperatures,  when  the  critical  current  is  lower  than  500  jiA,  the  I- Vs  were  found  to  . 
be  fit  quite  well  by  the  simple  resistively-shunted-junction  (RSJ)  model.  Upon 
cooling,  the  I-Vs  start  to  quickly  deviate  from  the  RSJ  model  and  show  an  excess 
current.  This  may  be  due  to  the  self-field  effect  when  the  width  of  the  junction 
becomes  much  larger  than  the  Josephson  penetration  depth  [17].  A  critical  current 
of  1  mA  corresponds  to  }c  =  1.2  x  105  A/cm2,  which  makes  the  Josephson  penetration 
depth  Aj  =  { ?i/2e^o(2A+i)/c}1/2  =  1  Jim.  This  value  is  smaller  than  the  junction 
width  of  4  jim.  However,  it  is  not  clear  whether  our  SNS  geometry  is  described  by 
an  overlap  geometry  or  an  in-line  geometry  as  discussed  by  Barone  and  Patemo  [18]. 
At  lower  temperatures,  when  the  critical  current  is  higher  than  ~5  mA,  the  I-V 
becomes  a  flux-creep-like  curve,  even  above  the  bulk  Tc  of  50  K  we  have  measured 
on  our  unpattemed  Ca-doped  YBCO  films.  We  found  that  the  crossover  to  flux- 
creep  behavior  depends  on  the  magnitude  of  the  critical  current  rather  than  the 
temperature. 
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This  junction  was  also  tested  at  77  K  in  a  magnetic  field  and  under 
microwave  radiation.  When  the  magnetic  field  is  applied  in  the  direction 
perpendicular  to  the  substrate,  we  observe  the  characteristic  diffraction  pattern  of 
IC(H)  with  a  central  peak  and  an  Ic  modulation  with  a  period  of  0.15  Gauss  as  shown 
in  Fig.  4.  We  do  not  understand  quantitatively  what  diffraction  pattern  should  be 
expected  for  the  edge  junction  geometry  because  of  the  overlap  of  the  second  YBCO 
film,  the  anisotropy  of  the  materials,  and  the  flux  focusing  effect  [19].  The  part  of  the 
second  YBCO  that  overlaps  the  junction  may  shield  the  initial  applied  magnetic 
field,  resulting  in  an  apparent  shift  of  the  magnetic  field  value  where  the  critical 
current  is  maximum.  Microwave  irradiation  modulates  the  critical  current,  and  all 
the  junctions  show  well-defined  Shapiro  steps  as  shown  in  Fig.  5. 

The  temperature  dependence  of  the  critical  current  for  junctions  with  200  A, 
400  A,  and  600  A  thick  barriers  is  displayed  in  Fig.  6(a).  The  two  sets  of  data  for  each 
barrier  thickness  represent  junctions  with  the  maximum  and  the  minimum  critical 
current  of  the  five  on  a  single  chip.  The  temperature  dependence  of  the  critical 
currents  shown  in  this  figure  is  completely  different  from  the  dependence  observed 
in  most  high-Tc  Josephson  junctions,  for  which  a  quasi-linear  dependence  over  a 
wide  temperature  range  is  usually  reported.  The  critical  currents  we  observe  are 
also  higher  than  the  critical  currents  reported  in  other  YBCO-based  junctions.  For 
instance,  the  critical  current  for  a  barrier  of  400  A  is  as  high  as  9  x  104  A/cm2  at  77  K, 
which  can  be  important  for  applications  requiring  a  high  }c. 

Figure  6(b)  displays  the  temperature  dependence  of  the  junction  resistance  Rn 
for  the  previous  junctions.  The  resistances  were  determined  by  the  slope  in  the 
linear  part  of  the  I-V  curve  at  a  bias  current  of  roughly  five  times  Ic.  The  samples 
with  barriers  of  200  A  and  400  A  show  a  linear  increase  of  Rn  with  the  temperature, 
but  the  data  for  the  600  A  sample  show  a  slight  structure  around  80  K.  This 
structure  is  highly  reproducible  over  several  junctions  on  different  chips  with  the 
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same  barrier  thickness.  We  also  note  that  the  resistance  does  not  scale  well  with  the 
thickness  of  the  barrier  at  high  temperatures  where  the  thinner  barrier  junctions 
have  higher  resistances  than  the  bulk  resistivity  would  predict.  The  dashed  lines 
are  the  expected  values  of  junction  resistances  for  200  A,  400  A  and  600  A  thick 
barriers  as  calculated  from  the  bulk  resistivity  in  Fig.  1  and  from  the  geometry.  At 
temperatures  near  75  K,  the  junction  resistance  (extrapolated  for  200  A)  appears  to 
scale  well  with  the  thickness.  Most,  interestingly,  it  is  obvious  that  the  junction 
resistance  decreases  faster  than  the  bulk  resistivity  would  predict.  For  all 
thicknesses,  the  resistance  is  too  high  at  high  temperatures,  correct  near  75  K,  and 
too  low  at  low  temperatures  compared  to  the  calculated  values.  As  far  as  we  know, 
this  behavior  has  never  been  reported  for  high-Tc  Josephson  junctions  which  have 
high  (10'8  Q-cm2)  interface  resistances.  However,  a  similar  behavior  has  been 
observed  in  the  Co-doped  YBCO  based  junctions  [10,  11,  13],  particularly  for  those 
with  low  doping  levels.  The  origin  of  this  temperature  dependence  of  Rn  is  still  not 
well  established,  even  though  Koren  et  al.  explained  this  behavior  by  the  flow  of 
Josephson  vortices  in  the  junctions  [11].  In  the  following  sections  we  will  offer 
alternative  interpretations  for  the  observed  temperature  dependence  of  Rn. 


IV.  Analysis  with  the  proximity  effect 


We  will  show  in  the  following  that  the  behavior  of  IC(T)  can  be  explained  if 
there  is  a  proximity  effect  in  the  Ca-doped  layer.  De  Gennes'  theory  for  a 
conventional  proximity  effect  predicts  in  the  dirty  limit  and  close  to  Tc  [20]  that 
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where  t  =  T/TC,  Iq  —  7tAq  /  AeRnkBTc,  d  is  the  thickness  of  the  barrier,  and  k_1  is  the 
decay  length.  In  the  case  of  a  BCS  gap,  2Ao  =  3.52  k^Tc,  and  for  Tc  =  89  K,  Io=  18.5 
mV/Rn.  In  the  absence  of  definite  knowledge  of  the  superconducting  gap  value,  we 
used  a  prefactor  that  will  generate  a  voltage  scale  close  to  a  presumed 
superconducting  gap  value  of  YBCO,  20  mV.  Even  though  the  relation  of  Eq.  (1)  was 
originally  derived  in  the  dirty  limit  near  Tc,  as  a  first  attempt  to  fit  the  general 
temperature  dependence  of  the  critical  current  over  the  entire  temperature  range, 
we  will  use  Eq.  (1)  as  a  generalized  form.  Since  the  electronic  mismatch  between 
YBCO  and  Ca-doped  YBCO  is  not  significant,  we  used  a  (1-f)2  dependence,  which 
assumes  continuous  boundary  conditions  of  the  wave  function  and  its  derivative 
[1].  However,  we  also  tried  to  fit  our  data  with  the  (1-f)  dependence  characteristic  of 
tunnel  junctions.  It  is  obvious  that  the  fitting  is  not  as  good  as  the  (1-f)2  case, 
especially  for  temperatures  close  to  Tc  where  the  (1-f)  or  (1-f)2  factor  is  dominant, 
since  in  this  range  the  fit  is  relatively  insensitive  to  the  length  dependence  [21].  This 
is  in  contrast  with  the  recent  analysis  of  our  data  [e.g.  that  in  Fig.  6(a)]  by  Kleinsasser 
and  Delin  [22],  who  used  a  (1-f)  dependence  by  assuming  rigid  boundary  conditions 
as  in  tunnel  junctions. 

When  the  barrier  material  has  a  superconducting  transition  at  Tn,  the  decay 
length  is  expressed  in  the  dirty  limit  by 


k  = 
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where  the  diffusion  constant  D  =  Vfl/3,  and  l  is  the  electronic  mean  free  path.  For 

the  clean  limit  the  decay  length  can  be  expressed  [1]  within  the  Ginzburg-Landau 
theory  by 


kt1  =  V2  x  0.74  ^0  (T/Tn  - 1)'1/2, 
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where  £0  =hvF/7tA  is  the  BCS  coherence  length.  We  have  used  a  temperature- 
independent  resistance  for  the  fitting,  since  the  origin  of  the  temperature-dependent 
resistance  is  not  clear  at  this  time.  We  will  discuss  this  issue  in  the  next  section  of 
this  paper.  As  in  the  Co-doped  YBCO  system,  the  fits  using  the  dirty  limit  [Eq.  (2)] 
and  the  clean  limit  [Eq.  (3)]  are  the  same.  This  is  not  surprising  since  it  is  easy  to 
show  that  Eq.  (2)  can  be  modified  to  K"1  (T/Tn-  l)'1/2  when  T  is  near  Tn,  where 
even  the  dirty  limit  formula  is  expected  to  follow  the  Ginzburg-Landau  relation.  In 
order  to  determine  in  which  limit  are  the  Josephson  junctions  with  Co-doped  YBCO 
barriers,  we  have  looked  at  the  variation  of  and  D  with  the  Co  doping  [13];  the 
data  are  consistent  with  the  barriers  with  low  Co  doping  (4%  and  7%)  being  in  the 
clean  limit,  and  the  barrier  can  be  described  by  the  dirty  limit  formalism  when  the 
Co  doping  is  14%.  However,  in  the  present  study  this  procedure  is  not  applicable 
because  we  have  studied  only  one  Ca  doping  level  owing  to  the  difficulty  of 
controlling  compensating  defects.  We  will  show  that  we  can  still  determine  the 
limit  of  the  barrier  by  comparing  the  Co-doped  and  the  Ca-doped  system. 

In  the  clean  limit  [Eq.(3)]  the  Ic  values  have  been  fitted  with  L  being  the 
only  parameter,  where  the  dimensionless  L  is  defined  by  Kd  =  L(t/tc)1/2  .  The 
fitting  lines  shown  in  Fig.  6(a)  have  been  calculated  with  Tn=  60  K  and  with  L  =  3.3, 
7.4,  and  12.9  (for  a  barrier  thickness  of  d  =  200  A,  400  A,  and  600  A)  which  correspond 
to  an  average  coherence  length  qo  =  52  A.  For  each  barrier  thickness,  L  is  an  average 
value  of  the  fitting  parameter  used  on  the  two  sets  of  data.  This  gives  a  decay  length 
as  high  as  100  A  at  77  K,  which  is  also  the  normal  coherence  length  we  estimated 
from  the  exponential  decrease  of  the  critical  current  with  the  thickness.  We  can 
compare  this  value  to  =  30  A  obtained  for  the  clean  4%  Co-doped  YBCO  which  has 

a  critical  temperature  Tc=  50  K,  similar  to  the  Ca-doped  YBCO.  From  these 
measurements  and  with  frtA,  the  ratio  between  the  Fermi  velocity  cy  (Ca) 
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and  Vf  (Co)  should  be  1.7  which  corresponds  to  a  ratio  of  3  (in  the  case  of  2D  system) 
or  5.2  (in  a  3D  system)  in  the  carrier  density  n  (if  we  assume  that  it  is  valid  to  use  the 
free  electron  model  to  estimate  n  from  Vf).  Carrington  et  al.  [23]  have  reported  Hall 
measurements  on  Co-doped  single  crystals  which  show  that  the  ratio  between  the 
carrier  density  of  YBCO  and  that  of  the  4%  Co-doped  YBCO  is  about  3.5.  For  the  Ca- 
doped  system,  the  carrier  density  is  sensitive  to  the  oxygen  concentration  [15],  but  if 
we  assume  that  the  only  reason  our  films  have  a  Tc  of  50-60  K  is  because  they  are 
overdoped,  then  we  can  roughly  estimate  the  ratio  between  the  carrier  density  in  the 
Ca-doped  YBCO  and  YBCO  by  using  an  empirical  formula  extracted  from 
thermopower  measurements  of  Obertelli  et  al.  [24]  and  confirmed  in  Ca-doped 
YBCO  thin  films  by  Kucera  et  al.  [15].  We  found  that  this  ratio  is  about  1.4,  implying 
that  the  ratio  n(Ca)/n(Co)  is  about  5,  which  is  close  to  the  value  we  have  extracted 
from  our  fitting. 

If  the  Ca-doped  YBCO  is  in  the  dirty  limit,  where 

2  T1/2 

1  h - ,  the  fitting  of  IC(T)  with  Eq.  (2)  and  Tn=  60  K  gives  L  = 

HT/Tn)J 

5.7, 12.7,  and  22  ,  corresponding  to  a  diffusion  constant  D  =  7  cm2/sec.  This  value  is  a 

factor  of  1.5  higher  than  the  D=  5  cm2/sec  found  for  the  dirty  14%  Co-doped  YBCO 
[13].  Since  D  =  Vfl/3  and  l  =  mv f  j (e2npn),  where  pn  is  the  resistivity  of  the  barrier,  it 

can  be  easily  shown  that 

D{Ca)  (n(Co)^/3pn(Co) 

D(Co )  U  (Ca)J  pn(Ca) 

where  the  Ca  index  is  for  the  Ca-doped  YBCO  and  Co  stands  for  the  14%  Co-doped 
YBCO.  The  Co-doped  YBCO  resistivity  pn(Co)  is  about  800  (iD-cm  at  100  K  compared 
to  60  jiO-cm  for  the  Ca-doped  YBCO.  From  the  data  of  Carrington  et  al.  [23]  we  can 
extrapolate  that  the  ratio  n(Ca)/n(Co)  is  about  13  and  then  the  ratio  D(Ca)/D(Co) 
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should  be  5.6,  which  is  higher  than  the  ratio  1.5  obtained  from  the  fitting.  The 
comparison  with  the  Co-doped  barrier  then  seems  to  indicate  that  our  Ca-doped 
barrier  is  indeed  in  the  clean  limit,  again  in  contrast  with  the  recent  analysis  of  our 
data  by  Kleinsasser  and  Delin  [22]  who  found  the  Ca-doped  YBCO  barrier  to  be  in  the 
dirty  limit.  Our  conclusion  is  also  confirmed  by  the  fact  that  the  resistivity  of  the  Ca- 
doped  YBCO  is  lower  than  YBCO,  which  is  known  to  be  in  the  clean  limit  [25]. 

Turning  now  to  the  anomalous  temperature  dependence  of  Rn  [Fig.  6(b)],  we 
note  that  a  somewhat  similar  observation  of  a  fast  decrease  of  junction  resistance 
has  been  reported  by  Clarke  et  al.  [26]  in  low-Tc  Pb-Ir-Pb  junctions  in  which  Ir  is  a 
superconductor  with  Tc  of  57  mK.  Following  the  argument  made  to  explain  the 
low-Tc  system,  we  suggest  that  one  explanation  of  the  decrease  of  Rn(T)  we  observe 
in  our  junctions  may  also  be  attributed  to  a  proximity  effect  in  the  Ca-doped  YBCO. 
According  to  the  theory  of  the  proximity  effect,  when  a  normal  metal  is  in  contact 
with  a  superconductor  through  a  very  clean  interface,  the  pair  potential  An  in  the 
metal  is  either  zero,  if  the  electron-electron  interaction  V n  =  0,  or  decreases 
exponentially  from  the  value  Ani  at  the  interface  if  Vn  >  0.  We  now  consider  a 
current  flowing  through  the  S-N  interface  and  an  electron  coming  from  the  N  side. 
Andreev  reflection,  the  conversion  of  normal  current  to  superconducting  current, 
occurs  at  the  point  where  An  is  equal  to  the  energy  of  this  electron.  If  Vn  =  0,  the 
Andreev  reflection  takes  place  at  the  interface,  and  in  this  case  the  resistance  of  the 
junction  is  expected  to  be  equal  to  the  bulk  resistance  of  the  N  layer.  This  was 
observed  for  Pb-Cu-Pb  junctions  [27].  On  the  other  hand,  when  the  barrier  itself  is  a 
superconductor  and  Vn  >  0,  the  Andreev  reflection  takes  place  inside  the  N  layer 
and  the  junction  resistance  should  be  less  than  the  bulk  N  layer  resistance.  For  this 
situation  Clarke  et  al.  developed  a  simple  model  for  a  conventional  proximity  effect 
and  found  that  the  resistance  increases  logarithmically  with  temperature,  which  was 
in  very  good  agreement  with  the  experimental  results  for  Pb-Ir-Pb  [26].  Flowever, 
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our  data  show  that  the  junction  resistance  varies  almost  linearly  with  the 
temperature.  This  difference  could  be  due  to  stronger  coupling  in  our  case,  where  Jc 
is  very  high,  in  contrast  to  the  case  of  the  Pb-Ir-Pb  junctions  where  the  coupling  is 
very  small  and  Jc  is  zero.  Although  Clarke  et  al.  were  able  to  quench  the  induced 
gap  in  the  barrier  at  a  few  (J.V,  we  did  not  observe  such  an  effect  up  to  a  couple  of 
hundred  |±V.  This  could  be  due  to  the  difference  in  the  interaction  potential  in  Ir 
(Tc=  57  mK)  and  Ca-doped  YBCO  (Tc=  50  K).  We  also  note  that  this  model  for  Rn 
cannot  explain  the  fact  [Fig.  6(b)]  that  for  all  thicknesses  of  the  barrier,  Rn  is  higher 
than  expected  for  temperatures  in  the  range  75  to  88  K. 

V.  Analysis  using  a  diffusion  model 

Another  possible  explanation  of  Rn(T)  is  that  there  is  diffusion  of  the  Ca  from 
the  barrier  layer  into  the  electrodes.  For  example.  Fig.  7(a)  shows  the  ideal. 30%  Ca 
doping  profile  uniformly  distributed  across  the  barrier  thickness  d  compared  to  a 
proposed  diffused  profile  in  which  the  Ca  distribution  is  a  Gaussian  with  a  full 
width  at  half  maximum  of  d.  Because  the  amount  of  Ca  is  the  same  in  both  cases, 
the  maximum  doping  for  the  Gaussian  profile  is  28.2%.  This  doping  profile  results 
in  a  Tc  distribution  throughout  the  barrier  as  shown  in  Fig.  7(b).  Note  that  in  this 
model  we  assume  a  uniform  distribution  along  the  width  of  the  junction,  i.e.  the 
diffusion  results  in  a  varying  Ca  composition  only  transverse  to  the  barrier.  This  Tc 
profile  has  the  following  consequences.  At  high  temperatures  T  close  to  To  if  the 
effective  barrier  thickness  is  defined  as  the  thickness  of  material  with  TC<T,  then  the 
barrier  is  thicker  than  the  nominal  deposited  Ca  doped  layer.  At  lower 
temperatures,  using  the  same  definitions,  the  barrier  is  thinner  than  the  deposited 
layer.  Hence  at  high  temperatures  the  junction  resistance  will  be  larger  than 
expected,  while  at  low  temperatures  the  resistance  will  be  smaller.  This  is  what  we 
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observe.  The  data  of  Fig.  6(b)  suggest  that  near  75  K,  the  barrier  thickness  is 
comparable  to  the  deposited  layer  thickness,  since  the  resistances  scale  appropriately 
with  deposited  thickness  at  that  temperature. 

In  this  diffusion  model,  we  still  rely  on  the  proximity  effect  for  the 
temperature  dependence  of  the  critical  current.  Since  the  effective  thickness  of  the 
barrier  depends  on  temperature,  and  Tc  is  not  uniform  transverse  to  the  barrier,  it  is 
difficult  to  calculate  the  exact  temperature  dependence  of  the  critical  currents. 
However,  we  do  not  expect  the  overall  shape  of  IC(T)  to  be  altered  significantly  from 
what  Eq.  (1)  predicts.  We  are  currently  in  the  process  of  investigating  the  interface 
diffusion  issue  with  transmission  electron  microscopy. 

We  can  further  extend  this  model  to  take  into  account  the  possibility  of 
doping  variations  in  the  lateral  direction  of  the  barrier.  We  discuss  this  "percolation 
model"  in  the  following  section. 

VI.  Analysis  using  a  percolation  model 

Although  the  IC(T)  data  of  the  Ca-doped  barrier  junctions  may  be  described  by 
conventional  proximity  effect  theory,  this  observation  is  only  one  of  many  relating 
to  our  devices  which  we  do  not  in  fact  completely  understand.  Examples  include 
the  temperature  dependence  of  the  resistivity  of  oxygen-deficient  and  doped  YBCO, 
the  role  of  anisotropy  in  the  edge  junction  devices,  a  detailed  understanding  of  the  I- 
V  characteristics  and  the  magnetic  field  diffraction  patterns,  the  nature  of  the 
interface,  non-reproducibility  of  Ic  and  Rn,  and  even  whether  YBCO  is  BCS-like  or 
not.  Thus,  a  high-Tc  proximity  effect  is  hardly  proved.  As  has  often  been  the  case 
in  the  high-Tc  arena,  we  turn  again  to  phenomenological  models  to  guide  our 
thinking,  particularly  for  improved  junction  devices.  We  are  mindful  of  the 
considerable  evidence  suggesting  that  oxygen-depleted  YBCO  is  inhomogeneous 
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with  respect  to  its  superconducting  properties.  Studies  have  shown  that  the 
superconductivity  throughout  microscopic  regions  of  carrier-depleted  thin-film 
YBCO  is  of  a  filamentary  nature  [9,  28].  Indeed,  this  is  the  reason  why  an 
elimination  of  excess  interface  resistance  in  high-Tc  junctions  was  originally  desired 
as  a  way  to  obtain  uniform  device  characteristics  [8]. 

Therefore,  a  very  important  question  regarding  the  SNS  junctions  is  the  role 
played  by  the  inhomogeneities  in  the  doped  YBCO  system.  We  have  already 
pointed  out  that  the  properties  of  Ca-doped  YBCO  are  particularly  sensitive  to  both 
the  Ca  concentration  and  the  oxygen  concentration.  If,  as  the  evidence  suggests,  the 
Ca-doping  or  the  oxygen  concentration  or  the  carrier  density  are  not  homogeneous 
inside  the  junction  barrier,  we  can  represent  this  barrier  by  a  set  of  grains  with 
different  Tc,  as  shown  in  Fig.  8  (these  grains  may  or  may  not  correspond  to  the  actual 
crystalline  grains  forming  the  thin  film).  As  we  can  see  from  the  particular  model 
in  Fig.  8,  when  three  adjacent  grains  that  traverse  the  barrier  become 
superconducting,  a  new  filamentary  superconducting  path  is  formed  which  is  a 
weak  link  (or  can  be  formed  by  weak  links)  and  which  gives  rise  to  the  critical 
current  of  the  junctions.  The  nature  of  the  weak  coupling  of  these  filaments  is 
puzzling,,  though  no  more  so  than  for  well-researched  but  still  ill-understood  grain 
boundary  junctions.  The  dependence  of  the  supercurrent  on  temperature  can  be 
predicted  within  this  percolation-like  model. 

We  derive  a  very  simple  relation  by  first  assuming  that  the  distribution  D(T) 
of  Tc  inside  the  Ca-doped  YBCO  barrier  is  linear  between  a  temperature  Tmin  and 
Tmax  with  the  simple  form 

®(T)  =  C(1  -  — — — )  (5) 

•‘max 
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which  is  represented  in  Fig.  9(a).  Assumption  of  a  different  distribution,  e.g. 
Gaussian,  does  not  significantly  affect  the  results  of  the  model.  The  value  of  D(T) 
shown  in  this  figure  represents  the  probability  of  finding  a  grain  with  a  critical 
temperature  of  T.  Then  the  probability  P(T)  to  have  superconducting  grains  at  this 
temperature  is  represented  by  the  shaded  area  between  T  and  T^x.  It  is  easy  to  show 
that  if 


P(T)=  \v{t)dt  and 
T 


max 

1=  Jo(r)df, 


this  probability  is 


(6) 


P(T)  = 


(Tmax  ~  T) 
(^max  —  ^min  ) 


2  ' 


(7) 


P(T)  is  shown  in  Fig.  9(b)  for  a  linear  distribution  of  Tc  between  Tm/n  and  Tmax  .  In 
order  to  have  a  superconducting  path  through  the  barrier,  n  aligned  grains  should 
be  superconducting,  where  n  =  d/ dgrain,  d  is  the  barrier  thickness,  and  dgrain  is  the 
grain  size.  In  this  percolation  model,  the  probability  to  have  n  grains  aligned  is 
Pn(T),  and  the  critical  current  of  the  junction  is  then 

IC(T)  =  Jc(T)APn(T),  (8) 

where  /C(T)  is  the  critical  current  density  and  A  =  8  x  10'9  cm2  is  the  area  of  the 
junction  in  our  edge  geometry.  We  will  assume  that  the  critical  current  density  is 
that  of  bulk  YBCO,  i.e.  it  decreases  linearly  with  temperature,  }c  =  Jo  (1  -T/Tc),  where 
Tc  =  90  K.  We  fix  Jo=  2  x  107  A/ cm2,  which  corresponds  to  a  reasonable  value  of  Jc  =  3 
x  106  A/ cm2  at  77  K.  The  results  of  the  fitting  with  the  percolation  model  are  shown 
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in  Fig.  10.  The  data  for  all  barrier  thicknesses  are  fit  nicely.  The  averages  of  Tmin  arid 
Tmax  for  each  Pair  of  curves  are  found  to  be  61.7  K  and  88.4  K  for  the  200  A  barrier, 
57.1  K  and  88.9  K  for  the  400  A  barrier,  and  54.7  K  and  82.2  K  for  the  600  A  barrier. 
From  these  values  we  can  estimate  that  Tc  in  the  Ca-doped  barrier  is  distributed 
from  Tmin  =  58  K  to  Tmax  =  86  K,  which  is  consistent  with  the  variation  in  the  value 
of  Tc  observed  in  the  Ca-doped  system  [15].  We  have  used  for  this  fitting  n  =  1,  2, 
and  3  for  the  200  A,  400  A,  and  600  A  barrier  thickness,  respectively,  which  means 
that  the  size  of  the  grains  in  the  model  is  roughly  200  A.  As  enough  grains  become 
superconducting  and  the  regions  of  filamentary  superconductivity  become  large,  the 
weak  coupling  will  be  lost  and  flux-creep  TV  characteristics  will  be  observed,  as  we 
have  mentioned  to  be  the  case. 

If  we  take  into  account  the  effect  of  a  different  Tc  for  each  individual 
superconducting  filament  on  the  temperature  dependence  of  the  critical  current,  the 
Ic  (T)  relation  will  be  described  by 


T 

Ic(T)  =  A  J 


dfn(rc) 

dTc 


(9) 


One  can  easily  see  that  Eq.  (9)  becomes  Eq.  (8)  if  one  ignores  the  dependence  of  Jo(  1- 
T /T c )  on  the  Tc  variation.  The  curves  obtained  by  fitting  with  Eq.  (9)  are  very 
similar  to  the  curves  presented  in  Fig.  10  and  calculated  with  Eq.  (8).  The  only 
difference  is  a  slight  change  in  the  two  fitting  parameters;  Tmfn  and  Tmax  increase  to 
72.7  K  and  89  K  for  the  200  A  barrier,  66.8  K  and  89.5  K  for  the  400  A  barrier,  62.2  K 
and  83.5  K  for  the  600  A  barrier. 

One  difference  between  these  devices  and  grain  boundaries  or  other  types  of 
high-Tc  junctions  is  the  greater  thickness  of  the  effective  "barrier"  giving  rise  to  the 
Josephson  behavior,  i.e.  the  entire  doped-YBtZO  interlayer  is  important,  rather  than 
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the  disordered  region  near  a  YBCO/YBCO  or  YBCO/  metal  interface.  The  regions  of 
filamentary  superconductivity  are  relatively  larger  in  the  present  case,  however,  and 
there  is  little  or  no  strain  to  prevent  complete  superconductivity  at  lower 
temperatures,  hence  the  sharp  rise  in  Ic  with  ( T-Tn ).  We  should  mention  that 
because  of  this  larger  length  scale,  this  simple  model  alone  may  not  explain  the 
Josephson  behavior.  A  proximity  effect  between  the  superconducting  grains  may 
still  play  a  role. 

We  can  also  explain  qualitatively  the  behavior  of  the  junction  resistance  with 
temperature  within  the  percolation  model,  in  terms  of  non-connected  but 
superconducting  grains.  When  the  temperature  is  decreased  from  90  K,  a  fraction 
of  the  grains  will  become  superconducting,  and  it  is  likely  that  the  current  will 
preferentially  flow  through  these  regions.  The  presence  of  these  superconducting 
grains  in  the  middle  of  the  current  path  will  therefore  effectively  result  in  lower 
resistance.  As  the  fraction  of  these  superconducting  regions  increases  with 
decreasing  temperature,  the  resistance  will  decrease  faster  than  the  bulk  resistance  of 
the  barrier.  This  is  consistent  with  the  observation  that  the  rapid  increase  of  critical 
currents  corresponds  to  the  points  where  the  resistances  extrapolate  to  zero  (see  Fig. 
6). 


VII.  Conclusions 

In  section  IH  we  reported  the  electrical  properties  of  SNS  Josephson  junctions 
based  on  YBCO  with  Yo.7Cao.3Ba2Cu307-§barriers.  We  found  that  the  critical  current 
of  these  junctions  has  a  temperature  dependence  which  is  not  the  usual  quasi-linear 
dependence  observed  in  most  high-Tc  junctions.  Furthermore,  the  junction 
resistances  Rn  decrease  faster  with  decreasing  temperature  than  the  resistance 
calculated  from  the  resistivity  of  Ca-doped  YBCO  thin  films.  In  sections  IV,  V,  and 
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VI  we  described  three  models  which  can  explain  the  data  presented  in  the  previous 
section.  We  showed  in  section  IV  that  the  IC(T )  curves  can  be  fitted  and  Rn(T) 
qualitatively  understood  within  the  de  Gennes  theory.  However,  the  Rrt(T)  can  be 
better  understood  qualitatively  if  one  assumes  diffusion  between  the  barrier  and  the 
electrodes  as  described  in  Section  V.  On  the  other  hand,  those  temperature 
dependencies  can  also  be  explained  by  a  very  simple  percolation  model  in  which  the 
junction  barrier  is  formed  by  grains  with  different  critical  temperatures,  as  discussed 
in  section  VI. 

At  this  point  it  is  not  clear  which  model  describes  the  real  behavior  of  our 
Josephson  junctions.  Certainly,  the  fast  increase  of  the  critical  current  alone  does 
not  seem  to  guarantee  the  proximity  effect,  especially  in  this  complicated  materials 
system.  We  also  note  that  the  proximity  theory  is  not  explicably  consistent  with 
some  of  our  other  observations,  and  the  diffusion  and  percolation  models  are 
phenomenological  only.  The  degree  of  homogeneity  of  the  barrier  may  determine 
which  mechanism  explains  the  junction  properties,  and  it  will  certainly  determine 
the  uniformity  and  transport  properties  of  the  devices.  It  is  possible  that  multiple 
mechanisms  are  involved. 

From  the  technological  point  of  view,  these  Josephson  junctions  with  Ca- 
doped  barriers  may  be  useful  for  applications  that  require  high  values  of  Jc. 
However,  a  serious  limitation  is  the  high  sensitivity  of  Tc  to  oxygen  concentration 
and/or  Ca  distribution,  which  could  affect  the  reproducibility  of  the  junctions. 
Other  limitations  include  an  excessive  Jc,  a  rapid  variation  of  Jc  with  T,  and  an 
insufficient  RnA  for  junctions  made  by  regular  optical  lithography  (i.e.,  junctions 
larger  than  2  pm  in  width). 
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FIGURE  CAPTIONS 

Figure  1.  Temperature  dependence  of  resistivity  for  a  1000  A  thick 
Yo.7Cao.3Ba2Cu307-5  thin  film  on  a  LaA103  substrate  measured  on  an  unpattemed 
film  of  1/4"  x  1/4"  size. 

Figure  2.  Schematic  cross-sectional  drawing  of  the  edge  junction  geometry.  Both 
of  the  YBCO  electrodes  are  2000  A  thick. 

Figure  3.  Current- voltage  curves  of  a  YBCO  /  400  A  Ca-doped  YBCO  /  YBCO  edge 

junction  at  several  temperatures.  The  width  of  the  junction  is  4  |im. 

Figure  4.  Critical  current  of  a  YBCO/  400  A  Ca-doped  YBCO/YBCO  edge  junction 
as  a  function  of  an  applied  DC  magnetic  field  at  77  K. 

Figure  5.  I-V  characteristic  under  3.5  GHz  microwave  radiation  and  at  77  K  for 
the  same  device  shown  in  Fig.  4. 

Figure  6.  (a)  Critical  current  and  (b)  junction  resistance  as  a  function  of 

temperature  for  junctions  with  Ca-doped  YBCO  barriers  of  thickness  200  A,  400  A, 
and  600  A.  In  (a)  the  solid  lines  are  the  fits  using  the  proximity  effect  theory 
described  in  the  text.  The  fitting  within  the  clean  limit  and  the  dirty  limit  are 
overlapping.  In  (b)  the  dashed  lines  represent  the  calculated  values  of  the  200  A,  400 
A,  and  600  A  thick  barrier  resistance  using  the  resistivity  of  the  Ca-doped  thin  film. 

Figure  7.  (a)  Gaussian  distribution  of  Ca  doping  throughout  the  barrier  region  of 

thickness  d  assumed  in  the  diffusion  model.  The  dashed  lines  indicate  the 
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nominally  ideal  Ca  profile,  (b)  The  consequent  distribution  of  Tc  in  the  barrier 
region  due  to  the  Ca  doping  profile.  A  Tc  of  75  K  obtains  at  the  thickness  d  of  the 
nominal  barrier. 

Figure  8.  Schematic  drawing  of  the  edge  junction.  In  the  percolation  model  the 
barrier  is  represented  by  a  set  of  grains  with  different  Tc.  The  shaded  grains 
represent  the  formation  of  a  filamentary  superconducting  path. 

Figure  9.  (a)  Linear  distribution  D(T)  of  critical  temperature  in  the  Ca-doped 

YBCO  barrier  used  in  the  percolation  model.  The  shaded  area  represents  the 
probability  to  have  grains  with  a  critical  temperature  at  or  above  T.  (b)  Probability 
P(T)  to  find  grains  which  are  in  the  superconducting  state  at  a  given  temperature  T. 

Figure  10.  Critical  current  as  a  function  of  the  temperature  for  junctions  with  Ca- 
doped  YBCO  barriers  of  thickness  200  A,  400  A,  and  600  A.  The  solid  lines  represent 
the  fit  of  Ic  with  the  percolation  model. 
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